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This thesis describes the conversion of a computer 
program from Fortran IV for the NOS 1.2 operating system of 
the CYBER 175 or CDC 6600 computer to Fortran IV compatible 
with the Naval Postgraduate School IBH 3033 system. The 
converted program, called FL027 , calculates the inviscid, 
three-dimensional transonic potential flow over wings or 
wing-body combinations. The data input to FL027 is 
extensive; therefore, an interactive program was developed 
to aid the user in building the required input data file. 



4 



TABLE OF CONTENTS 



I. I NTEODUCT ION 9 

A. BACKGROUND 9 

B. 7ISC0CS/IN VISCID SYSTEM OF PROGRAMS 10 

II. POTENTIAL FLOW PROGRAM FL 027 14 

A. RE-PROGRAMING 14 

B. PROGRAM DESCRIPTION 19 

1. Program Input 22 

2. Program Output 32 

III. INTERACTIVE INPUT PROGRAM FL027IN 35 

IV. FL027 EATCH SYSTEM EXECUTION 38 

V. PROGRAM TEST RESULTS 40 

A. ACCEPTANCE TEST DATA 40 

B. COMPARISON WITH OTHER PROGRAMS 41 

C. AE-4501 CLASS PROJECT 41 

APPENDIX A 44 

APPENDIX E 46 

APPENDIX C ' 47 

APPENDIX D ’ 48 

APPENDIX E 70 

APPENDIX F 132 

LIST OF REFERENCES 151 

INITIAL DISTRIBUTION LIST 152 



5 



LIST OF T ABL ES 



I. Viscous/Inviscid Wing System of Programs .... 11 

II. C DC Magnetic Tape Files 15 

III. FL027 Re-Programing Changes 16 

IV. Initialized Input Variables 18 



6 



LIST OF FIGURES 



1.1. Viscous/In viscid Interaction Procedure 12 

2. 1. Dihedral Angle 27 

2. 2. Wing Defining Geometry 28 

2.3. Section Defining Geometry 31 

2.4. Versatec Plo*. of Cp vs. X/C 34 



5. 1. Program Calculated and Wind Tunnel Data .... 42 



7 



AC KNCWLEDG EMENT 



I would like to acknowledge the most valuable assistance 
of the W. R. Church computer center personnel in the 
successful completion of this thesis. In particular, the 
guidance of Mr Dennis Mar was instrumental in directing my 
initial efforts. I am grateful to Dr Robert Zucker for 
providing the constant attention and encouragement which 
kept me working at an even pace. In addition, I am deeply 
appreciative cf the love and support provided by my wife 
Tricia and our two children. Most of all I thank God for 
the strength, determination and peace which has sustained me 
and brought this project to a successful completion. 



8 



I. 



INTROD OCTION 



In the Aeronautical Engineering curriculum graduate 
level aerodynamics course, AE-4501, the sruder.ts are exposed 
to two computer programs. One of these, prepared by the 
Douqlas Aircraft Company, analyzes the potential flow around 
three-dimensional wings but is limited to incompressible 
flow THef. 1j. The other program, prepared by Cebeci, 
calculates the friction drag for two dimensional 
incompressible flow over airfoils [Ref. 2]. A serious 
defect of these programs is that they are not 
state-of-the-art computer programs. The Douglas program 
does not consider the effects of compressibility and the 
boundary layer program, in addition to being restricted to 
incompressible flow, does not predict the laminar to 
turbulent transition point. 

A. EACKGRCOND 

In 1980 the Department of Aeronautics at the Naval 
Postgraduate School acquired an intricate computer program 
recently developed by the Boeing Commercial Airplane 
ComDany. This state-of-the-art program calculates 
three-dimensional transonic flow over wings and bodies in 
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both the cuter-inviscid flow region governed by the 
transonic potential eguaticn and the thin layer in which the 
first order, compressible boundary layer equations are 
assumed tc te valid. 

The Boeina program as received was designed to be 
executed on a CDC 6600 or a CYB2R 175 computer and was 
written using CDC FORTRAN IV extended language. This thesis 
therefore was primarily concerned with the conversion of the 
oroqram to FORTRAN IV extended compatible with the Naval 
Postgraduate School's (NPS) IBM 3033 system. The large 
modular orcgram was divided so that the potential flow 
analysis portion could be run separately. Simplified 
instructions for use of the program were also prepared. 

3. VISCOOS/INVISCID ‘SYSTEM OF PROGRAMS 

The Viscous/Inviscid Wing System (VIWS) of programs 
calculates three-dimensional transonic flow over wings and 
wing bodv combinations including details of the laminar or 
turbulent flot. in the three-dimensional viscous boundary 
layer. The flew field is calculated in two overlapping 
regions: an cuter inviscid flow region governed by the 

transonic potential eguation, and a thin boundary layer in 
which the first order, three-dimensional, compressible 
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boundary layer equations are assumed to hold and in which 
the effects of surface heat and mass transfer can be 
computed. A list of the VIWS of programs is presented in 
Table I. 



TABLE I 

Viscous/In viscid Ming System of Programs 



Program Name 


Description 


F 1027 


Jamescn-Caughey inviscid, 
transonic wing code 


A 4 1 1 1N 


Reads geometry & velocity 
data, constructs coordinate 
system 


V WIN 


Potential flew boundary 
layer interface 


A411 AC1 


Three-dimensional boundary 
layer program 


INTERP 


Boundary layer potential 
flow interface 


A41 1E1 
A 4 1 IF 2 
A41 1 ES 


Graphics display programs 



The basic sequence of calculations used by the VIMS to 
obtain matched viscous and inviscid solutions consists of an 
iterative loop in which the inviscid outer flow analysis and 
the boundary layer analysis are performed sequentially. The 
iterative seguence is continued until either convergence 
(satisfactory matching) is achieved, or the maximum number 
of iterations specified by the user has been performed. The 
VIMS programing sequence is shown schematically in Fig. 1.1. 
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figure 1.1. Viscous/ In vise id Interaction Procedure 



The potential flow is 
the first iteration, 
of the boundary layer 



calculated for the bare wing during 
In subsequent iterations, the effect 
flew on the outer inviscid flow is 
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felt as a modification to the wing shape through the 
addition of the boundary layer displacement thickness. 
Convergence is recognized, and the iterations are stopped, 
when the maximum change between the new and old displacement 
thickness, expressed as a fraction of the maximum 
displacement thickness, is less than the convergence 
tolerance chcsen by the user. 

The VIWS utilizes the Jameso n-Caughey transonic inviscid 
wing program FL027, to carry cut the potential flow 
analysis. The boundary layer analysis is performed by a 
finite difference boundary layer prediction program 
developed by the Boeing Ccmmerciai Airplane Company. The 
basic theorv behind the bcundary layer program is contained 
in [ Hef . 3]. A detailed description of the VIWS of programs 
(excluding the potential flew program FL027) is contained 
in [ Sef . 41. A basic guide to the use of the VIMS of 
programs is contained in [Bef. 5]. 
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II. POTENTIAL FLOW PROGRAM FLQ27 



Eecause of the extensive length and number of program 
modules in the VIWS , the Potential Flow Program, FL027, was 
singled out for conversion. It was anticipated that FL027 
would be run separately at first and recombined with the 
other urogram modules at seme later date when zhese modules 
were themselves translated for execution on the IBM 3033 
computer. 

A. RE-PRCGBAMING 

The Potential Flow Program, hereafter called FL027, was 
received on magnetic tape and loaded into the IBM 3033 mass 
storage system using the Job Control Language (JCL) routines 
presehred in Appendix A. The magnetic tape contained twenty 
(20) tonal files in which the format was 9 track, 1600 CPI, 
unlatled. The card image format for the sixteen (16) 
program files is 80 characters per record and the four (4) 
output files contain 150 characters per record. The program 
and out out files on the original CDC nape are listed in 
Table II. 

The FL027 program was converted to FORTRAN 17 extended 
suitable for execution on the IBM computer using the NPS CDC 
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TABLE II 



C DC Magnetic Tape Files 



File/Records 


Name 


Description 


1 /2 356 


FLO 27 


Potential Flow Program 


2 /3 1 94 


A 4 1 1IN 


Reads geometry S velocity data 
const ruczs coordinate system 


3 / 378 


VWIN 


Potential flow boundary-layer 
inter face 


4 /6611 


A 4 11A01 


Three-dimensional boundary 
layer program 


5 /I 9 77 


INTERP 


Bound ary -layer potential flow 
inter face 


6 / 668 


A 41 IPS 


Streamline plots 


7 / 211 


A 4 1 IP 1 


One-dimensional plots 


8 / 586 


A 4 1 IP 2 


Contour plots 


9/70 


COUPLE 


Procedure files 


10 / 158 


ITER 


11 / 7 


DATAIN 


12 / 78 


FINAL 


12 / 434 


BOEB1 


Boeing McLean computer program 


14 / 36 


CONTPLT 


Contour plots 


15 / 17 


CORDPLT 


One-dimensional plots 


16 / 40 


STREPLT 


Strea mwise plots 


17 


OUTF27 


Output from FL027 


18 


OUTIFC 


Output from VWIN 


19 


0UT411L 


Output from boundary-layer, 
lower surface 


20 


OUT 41 1 U 


Output from boundary-layer, 
upper surface 



to IEM conversion guide [Bef- 6 ]. The f irsz step taken 
consisted cf program compilation using the WATFIV compiler 
with its extended error messages. The listing which was 
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produced flagged all araas of the program which required 
revision. Program changes were accomplished utilizing this 
WATFIV listing. Some of the mere general and repetitive 
chanqes are listed in Table III. 

TABLE III 

FL027 Re-Programing Changes 



CDC Code 


IBM Code Change 


Variables: FREAD, FREAF, 
FNRIT, FWR IF , IR E AD, 
IREAF, IWRIT, IWRIF 


Eliminated from program 


WRITE (I NR IT, 6 00) 


WRITE (6 ,600) 


R EAD ( IR EAF , 500) 


READ (5,500) 


READ 7, WRITE 7 or 
REWIND 7 


Changed to READ 14, 
WRITE 14 or REWIND 14 


Call SECOND (T) 


Step eliminated 


Call SSWITCH (1 ,ISTCP) 


Call SLITET ( 1, ISTOP) 


Delimiter of form * 


Replaced by * 


Comment cards with" * 


Replaced by C 


LEVEL statement 


Step eliminated 


If (UNIT (N ) .GT.O. ) GO TC 


All of this type eliminated 



The most 
Buffer IN or 
proqram to 
into and out 
memory size 
t-r :e cccurre 
subroutines. 



difficult change to make occurred with the CDC 
Buffer 00T statements which were used in rhe 
transfer portions of a three-dimensional array 
of main memory. The Buffer routines reduce the 
required to execute the program. This statement 
d in the main program and several of the 
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The change required to translate this statement is presented 
below with the CDC code preceeding the IBM FORTRAN. 

EOFFER OUT (N3,1) (G ( 1 , 1 , 1) , G ( M X , MY , 1 ) ) changed to 
WRITE (N3) ((G(I, J,1) ,1 = 1, MX) ,J*1,BY) and 

EOFFER IN (N 1 ,1 ) (G ( 1 , 1 ,M) , G (MX , MY , M) ) changed to 
READ ( N 1 , ERR= ) ( (G (I , J,M) ,1 =1 ,MX) , J= 1 ,MY) 

The variable ERR was assigned the GO TO statement number 
of the ONIT statement immediately following the BUFFER IN 
line of code. As an example, if the UNIT statement 
folic vino the EUFFER IN cede was - If (UNIT (N 1 ) . GT. 0. ) GO TO 
151, then the number 151 was assigned to variable ERR 
follcwinq the equal sign. All CDC UNIT statements were 
eliminated from the FL027 source code par Table II. 

In addition to the pregram changes required to run FL027 
on the IBM computer, several lines of code were added to 
modify the output format to a more usable form. A 
subroutine, VERTEC, which calls the Versatec plotter was 
also added to enhance program usefulness. This plotting 
routine is user controlled through an input variable and is 
ex lained in the next section. The modified FL027 program 
source code is presented in Appendix E. 
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Tc facilitate program daxa entry several inpur. variables 
which had recommended values were initialized to these 
values within the Main program and the subroutine GEOM. The 
initialized input variables and xheir values are presented 
in Table IV. 



TABLE IV 

Initialized Input Variables 



AREA 


VARIABLE NAME 


INITIALIZED VALUE 




XSCAL 


o 

• 

o 




PSCAL 


o 

• 

o 




FCONT 


o 

• 

o 




P20 


0.7 


MAIN 


P30 


1 .0 


Pr gm . 


FSMCO 


o 

• 

o 




PTMAP 


o 

• 

o 




BLCP 


o 

• 

o 




WEIG 


1 .0 




PTC K 


o 

• 

o 




FIX 


0.0 




YSYM 


o 

• 

o 




FNB 


2.0 




?X 


o 

• 

o 


Sufcrt. 


PZ 


0 .0 


GEOM 


TRL 


o 

• 

o 




SLT 


o 

• 

o 




XSING 


o 

• 

o 




YSING 


o 

• 

o 
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A complete description of each input variable in Table 
IV can be found on pages 19 through 23 of [Ref. 5]. 

B. PROGRAM DESCRIPTION 

The FLC27 program is a computer code written to analyze 
the transonic flow over a wing alcne or a wing on a 
cylindrical fuselage. It uses a finite-volume formulation 
to solve the exact potential flow equation in conservative 
form. In the development of the equations, the basic 
assumptions are; steady flow, no heat or work transfer, 
isentropic flow, irrotational flow, no body forces and a 
oerfect qas. The velocity vector in cartesian coordinates 
is 



where u, v and w are the velocity components. The 
continuity aquation, assuming steady flow, is 



Next a velocity potential is introduced such that the 
velocity components are calculated as the gradient of this 
oo tential . 



T /y /% 

= ui + vg + wk 



( 2 . 1 ) 




( 2 . 2 ) 



u 




(2.3) 
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With the introduction of the velocity potential, 
continuity equation 2.2 becomes 



the 



^t«^> + - 0 < 2 -“> 

Assuming no heat or work transfer, the energy equation can 
be written as 



♦ <r * ” » 2 1 


= T 


-Jk 

♦ 

1 

—A 

BE 

1 

1 


2 


«Q 


2 



The flow is assumed to be uniform in the far field. On the 
surface of the body, the normal velocity component is zero. 
The velocities and densities of the near field are 
normalized using the free stream velocity and density, thus 
= 1 and p. = 1 . Using the assumptions than the flow is 
isentropic and a perfect gas, the energy equation 2.5 can be 
shown to be 



P 



( 7 - 1 ) 2 f 

+ 2 (rv 



7-1 



( 2 . 6 ) 



With equations 2.5 and 2.6 there are only two unknowns, 
(p and p . They can be solved, subject tc the boundary 
condition of flew tangency, using a finite volume technique. 
The basic numerical scheme for the solution is the 
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construction cf a mesh ficm small volume elements (cubes) 

/ 

which are packed around the wing or wing body configuration. 
The cubes in the computational domain are separately mapped 
to distorted cubes in the physical domain by independent 
transforma tions from local coordinates X, Y and Z to 
Cartesian coordinates x, y and z. The mesh points are the 
vertices (corners) of the mapped cubes. The velocity 
potential and density are calculated an each vertex in the 
mesh. The pressure distribution can then be calculated from 




In the event that the local flow velocity becomes 
supersonic and shocks occur, these are handled in the usual 
manner bv insuring that: 

1) The tangential-' velocity components are equal on each 
side cf the shock. 

2) Continuity is maintained by keeping the product of 

'P 0^ constant across the shock (where 0^ is the normal 
velocity component) . 

3) Discontinuous expansions (corresponding to an 
"expansion shock") are excluded from the flow field. 
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The assumption cf isentropic flow along with the existence 
of shocks cresents a contradiction which can only be 
resolved by limiting the flow to very weak shocks for which 
entropy and vcrticity generation may be ignored. Thus, 
solutions will be valid only for subsonic free stream 
velocities. 

The main three-dimensional array containing the 
potential function data is stored on disk, and special 
unformated incut/output statements are used to bring planes 
of data into central computer memory and to store updated 
planes of data back on the disk. In the construction of the 
computational coordinate system, a Joukowski transformation 
is used to transform the cylindrical fuselage to a vertical 
slit and then a sheared parabolic transformation is used in 
planes containing the airfoil sections. A detailed 
mathematical formulation of the potential flow analysis is 
contained in [Hef. 7]. 

1 . Program Input 

The incut to FL027 consists of variables which are 
read with an 8F10.6 FORMAT. Each input card has a title 
card which precedes it. This title card contains the input 
variable name and effectively labels the input data for easy 
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reference. The title for each input data is placed in the 
same column as the input data it labels. The title cards 
are read with a 20A4 FORMAT. All numerical input values are 
real numbers. The following data deck, listed card by card, 
is the minimum input data required for a simple wing 
analysis. Bach "card" can be interpreted as one line of 
data on your terminal. A complete sample data sen is 
presented in Appendix C. 

CARD 1 The Run Title (64 characters maximum) 

CARD 2 Title card for the input variables 
FNX, FNY , FNZ , FMESH and FPLOT 

CARD 3 

Cols. 1-10 FNX - Number cf computational cells in 
the chordwise direction for the initial mesh. 

MAX = 160/2**n, where n = FMESH -1. (See Cols. 31-40 
for FMESH) 

Cols. 11-20 FNY - Number of computational cells in 
the normal direction frcm the airfoil surface for 
the initial mesh. 

MAX = 16/2**n, where n = FMESH -1. 

Cols. 21-30 FNZ - Number of computational calls in 
the spanwise direction for the initial mesh. 

MAX = 32/2**n , where n = FMESH -1. 
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Cols. 31-40 FMESH - Determines the number of times a 



CARD 4 



CARD 5 



program generated computational mesh is refined. 
Enter only 1 . 0, 2.0 or 3.0 for coarse, medium or 

fine mesh. If 3.0 is selected the program will 
calculate flow over the wing for the coarse mesh 
then half the mesh size (medium) , recalculate, then 
half the mesh again (fine) and do a final potential 
flew calculation. Output parameters are printed for 
each mesh size fer which calculations were 
performed. 

Cols. 41-50 FELOT - Output flag 

0.0 = Normal output without printer-plot of Cp 

1.0 = Normal output with printer-plot of Cp 

at each computational mesh point for 
each wing section. 

2.0 = Normal output with Versatec plots of 

Cp versus X/C for each wing section of 
the final mesh. 

Title card for the input variables 
FIT, COVO and P10 

-M Cne card for each computational mesh. Total 
number of cards egual to M = FMSSH. 
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CARD 6 



CARD 7 



Cols. 1-10 FIT - A parameter which fixes the 
maximum number of iterations the program will use to 
converge the velocity potential to a specified 
tolerance (COVO) . This parameter must be repeated 
for each mesh refinement. 

Cols. 11-20 COVO - Velocity potential convergence 
criteria. This input variable is also entered for 
each selected mesh. A value of 0.000001 is 
recommended . 

Cols. 21-30 P10 - This parameter determines the 

subsonic point relaxation factor for the specified 
mesh size. A value of less than 2.0 must be entered 
for each designated mesh. Recommended values are: 
1.6 for coarse, 1.3 for medium and 1.2 for the fine 
mesh. 

Title card for the input variables 
FMACH, YA, AL and CD0 

Cols. 1-10 FMACH - Free stream Mach number 
Cols. 11-20 YA - Yaw angle in degrees 
Cols. 21-30 AL - Angle of attack in degrees 
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Cols. 31-40 CDO - Drag coefficient due tc skin 
friction. Onless known, an estimated value of 0.01 
is rsccmmen ded. 

CARD 8 Title card for the input variables 
ZSYK, FNS , SWEEP, EIHED and FUS 

CARD 9 

Cols. 1-10 ZSYfJ - The wing planform symmetry 
trigger. 

0.0 = Yawed wing, has no spanwise symmetry 
1.0 = Swept wing, has spanwise symmetry 
Cols. 11-20 FNS - This input variable cells the 
orcgram the total cumber of wing sections you have 
selected to define the wing half span. The number 
must be at least three (3) but not more than eleven 
( 1 1) sections . 

Cols. 21-30 SWSEE - Leading edge sweep angle in 
degrees. 

Cols. 31-40 DIHED - Dihedral angle in degrees. See 

Fig. 2. 1. 

Cols. 41-50 FDS - Input the fuselage radius. Enter 
0.0 for a wing-alone case. 
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Data input cards from 1 0 through 15 are used for 

defining wing planforms and section geometries. For the 

€> 

first wing section, all data cards from card 10 through card 
15 must be used. For the second and subsequent sections - 
there is an option for skipping the wing section defining 
daxa (cards 12 through 15) and copying -he data from that of 
the previous section. This option is controlled by the input 
variable FSEC. If this option is not used, data cards from 
10 through 15 must be repeated for each wing defining 
section. The number of wing sections which are defined is 
input with the variable F NS. Remember, up to 11 sections may 
be defined, and a minimum of 3 sections must be defined. 
Ml wing planform and section defining geometries must be in 
consistent units. Wing planform and section defining 
quantities are presented in Fig. 2.2. 



27 




Figure 2.2. Sing Defining Geometry 

CARD 10 Title card for the input variables 

ZS, XI, YL, CHORD, THICK, AT and FS EC 

CARD 11 

Cols. 1-10 ZS - The section spanwise coordinate 



(Srart at the centerline and work 
cut bo ar d) 



Cels. 11-20 


XL - 


Sectio n 


leading 


edge 


X coordinate 


Cels. 21-30 


YL - 


Section 


leading 


edge 


Y coordinate 


Cols. 31-40 


CHORD 


- Sect; 


.on chord 


length 
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Cels. 41-50 THICK - The thickness scaling factor 
can be used to scale all Y coordinates of the wing 
section. Thus percent thickness and camber are 
increased (or decreased) accordingly. Use 1.0 if 
no scaling is desired. 

Cols. 51-60 AT - The twist angle of each section 
(geometric twist) measured from the X axis to the 
chord line. A positive twist angle reduces the 
section angle of attack and gives "washout". Use 
0.0 fer no twist. 

Cols. 61-70 FS EC - This is a flag which determines 
whether or not the program reads wing section 
defining geometry from a previous wing section or 
from new defining geometry. For the first section 
defined you must set FS EC to 1.0. Following the 
first section, if you define new section geometry 
then use FSEC = 1.0. If you want the program to 
read the section geometry defined from the previous 
section then set FSEC = 0.0. 

CARD 12 Title card for the input variable 
FN 
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CARD 13 



Cels. 1-10 FN - This variable contains -he number 



cf points which define the upper and lower surface 



of the section. A maximum of 161 points may be 



used. 

CARD 14 Title card for the input variables 
XF (I) and YP (I) 

CARDS 15-1 to 15-N Total number of cards equals N, 
where N = integer part of (FN+2)/3. 

The X and Y coordinates at each point are entered 
in Dairs, three points to a card. (See Appendix C 
for samcle input) 

Cels. 1-10 XP (I ) - X coordinate of the wing 

section point 

11-20 YP ( I ) - Y coordinate of the wing 
section point 



21-30 



defining X coordinate for next 



poi nt 



31-40 



defining Y coordinate for next 



poi nt 



4 1-50 



defining X coordinate for 



following point 
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5 1-6 0 



defining Y coordinate for 
following point 



The X and Y coordinates of the wing section defining 
points must be entered starting with the upper surface 
trailing edge pcint and proceeding along the upper surface 
ro the leading edge, and returning along rhe lower surface 
to the lower surface trailing edge point. It is very 
im ortant to define the section leading edge with a large 
number of closely spaced points. Suggest at least 0.05 
soacing or less between X coordinate values from 0.1 X/C to 
the leading edqe, X/C = 0.0. Each X and Y coordinate point 
is normalized using the chord length for that section. 
Section defining geometries are illustrated in Fig'. 2.3. 




Figure 2.3. Section Defining Geometry 
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CARD 16 Title card containing the words in Cols. 1-80 
END OF CALCULATIONS 

CARD 17 Title card for the input variable 
FNX 

CARD 18 

Cols. 1-10 FNX - This variable indicates the end 
of a set of calculations and must be set equal to 
0.0. Its purpose is to indicate that the program 
has run to completion. 

2. Program Out put 

Output from the FLC27 program varies with the value 
of the input variable FPLCT. When FPLOT is set equal to 0.0 
a normal output is produced. This normal output contains 
(in order of cccurance) : refined input geometry data 

including trailing edge slope and angle calculations ; 
iterative solution of the potential flow mesh; section 
characteristics and wing characteristics. The iterative 
solution, section and wing characteristic data are repeated 
for each mesh refinement requested. Thus, if the input 
variable FMESH is set equal to 3 , these data are calculated 
and output three times. The last data in the normal output 
consists of the non-dimension ali zed chord (X/C) and pressure 
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coefficient (Cp) data at each computational mesh poinc fcr 
each winq section calculated during the final mesh. A 
sample of the normal output data is presented in Appendix D 
and recresents the output data frcm Appendix C input data. 

If variable FPLOT is set equal to 1.0, the output 
data is increased considerably. This output contains the 
normal output plus a line printer-plot of the«pressure 
coefficient at each computational mesh point for each wing 
section. The line printer-plot is produced for each wing 
section of each mesh refinement. The length of the output 
data with FPLOT set equal to 1.0 can approach 6000 records 
defending on the number of mesh refinements requested. 
These clots are of questionable value and, therefore, an 
alternate plotting program was developed. 

When the variable FPLOT is set equal to 2.0, the 
normal output data is produced plus a Versatec plotting 
subroutine (VERTEC) is called. The subroutine outputs, via 
the Versatec plotter, plots of Cp versus X/C for each wing 
section of the final mesh calculations. This routine is 
simply putxinc into plot form the Cp and X/C numerical data 
contained in the normal output. A sample of the Versatec 
plot is oresented in Fig. 2.4. 
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SECTION CP DATA 
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Figure 2. 4, Versatec Plot of Cp vs. X/C 
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III. INTERACTIVE INPUT PROGRAM FLQ27IN 



The input data file required for the FL027 program is 
extensive. Errors in input data FORMAT will cause program 
errors at execution time. In order to eliminate these 
errors and reduce the input data workload, a computer 
terminal interactive program was written. This interactive 
proqram, called FL027IN, is a user-friendly way of creating 
an input data file fcr the potential flow wing analysis 
proqram FL027. The FL027IN program source code is presented 
in Appendix F. 

The interactive program, FL027IN, when executed ask's 
questions cf the user in order tc construct and write to the 
user's "A" disk the required FL027 input data file. The 
following presents the step-by-step procedure for executing 
the interactive program FL027IN. 

STEP #1 Log on to any IBM 3033 interactive terminal 

with your user number and password. 

STEP #2 Once logged on and in the CMS operation mode 

type : 

CP LIBK 0247P 191 120 RR then hit ESTER 

STEP #3 The word PASSWORD will appear. Type and ENTER 

AERO 
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STEP #4 Type and ENTER 

acc 120 D 

STEP #5 Type and ENTER 

LOAD PL027IN (START 

The screen will display the header for the interactive 
program. Answer each question presented. At the end of 
each question in parenthesis is the input data variable 
associated with that question and whether the input 
tarameter is a real number (R) or an integer (I). Example: 
==> Enter the free stream Mach number (FMACH) : (R) . FMACH 

is the input data variable for the question. As you proceed 
through the FLC27IN program, opportunities to review and 
change input data will be presented. Should it become 
necessary to change your input data after completing the 
FL027IN program, you can simply XEDIT the created data file. 

The FL027IN program also incorporates a library which 
contains the wing-section defining data for a number of 
current winq shapes. A copy of this library is presented in 
Ao'endix B. This feature will be displayed during program 
execution by the use of a menu from which the user can 
select a pre-defined wing section or define his own. 
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Upon completion of user inputs to the interactive 
oroaram three additional data lines are automatically 
written to the bottom of the input file. They are: 

END OF CALCULATION 
FNX 
0 . 0 

In addition. Job Control Language (JCL) cards are 
written to the top and bottom of the file. All JCL cards 
start with a // format. After FL027IN has run to completion 
tv n e and enter RELEASE 191 to release the aero disk which 
was linked while executing the FL027IN program. The created 
data file is written to the user’s "A" disk with <filename> 
<filetvoe> of FL027 DATAIN. Additional changes can be made 
simply by entering the XEDIT mode and editing the file. 
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IV. FLO 27 BATCH SY STEM EXECUTION 



The ootential flow program FL027 can be executed after 
the inout data file has been created. The batch processor 
is required for FL027 execution because of the extensive CPU 
time needed to run the program. While in the XEDIT mode, a 
standard JCB card must be added to the top of the FL027 
DATAIN file prior to submission for job execution. The JOB 
card has the form: 

//jotname JOB (nnnn, pppp) , 'ident ' ,CLASS=J 
jobname = may contain up to 8 alphanumeric characters, 
the first of which must be alphabetic, 
nnnn = your user number 

pppp = project number, assigned by professor 
•ident' = contains the user's own identification 

information. A maximum of 20 characters may be 
contained within the single quotation marks. 

After adding the JOB card to your data file, you are 
ready to execute the program. Type SUBMIT FL027 DATAIN and 
press ENTER. Eatch runs are normally net worth waiting for. 
To inquire about the status of the job, enter INQ and the 
job name used on the JOB card or "logoff". If the system is 
busy and the naximum mesh size was selected, it may be 
several hours before your job is run. 
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When the jot is run the output will be spooled to the 
batch printer located next to the VM printer in the main 
computer building. The title at the top of the printout for 
batch lobs is the name entered on the JOB card. If it is 
desired to have the program output data spooled directly to 
the terminal, it will be necessary to add one additional JCL 
card to the input data set. This card must be placed 
immediately following the JOB card and has the form: 

//♦MAIN 0RG=NPGVM1. nnnnP 
where nnnn = ycur user number 

Inserting this card in the input data will cause all 
program output to be spooled to the user's virtual reader 
where it may he looked at, printed or transferred to his "A" 
disk. To enquire as to whether information is in the reader 
simply type RDR and hit enter, then follow the instructions 
on the screen. 
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V. 



PROGFAM TEST RESULTS 



The FLC27 program was tested in three stages; (1) during 
the reprograming phase for conversion completeness, (2) 
after sucessful conversion with suitable wing data for 
program accuracy and (3) during an AE-4501 class project. 

A. ACCEPTANCE TEST DATA 

To test and ensure that the FL027 program was converted 
to I EM compatible Fortran without error, an acceptance test 
data set was used. The acceptance test input and output 
data was supplied with the original CDC program source code. 
After conversion of the FI027 program to Fortran suitable 
for the NPS IEH system, the acceptance test input data were 
run and the output results compared to the output generated 
by the CDC system. 

Eoth output data sets were numerically exact when the 
FL027 program was run in double precision on the IBM system. 
If the program was run ir single precision, the numerical 
output values were exact tc the third decimal place. The 
difference in single precision accuracy occurs because the 
CDC system uses a 64 bit word length while the IBM system 
word length in single precision is only 32 bits. It was 
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decided that the IBM single precision accuracy was 
satisfactory . 

B. COMPARISON WITH OTHER PROGRAMS 

The FLC27 program was also tesred for accuracy by using 
the wing olanfcrm and section data from a NACA 572 wing. 
The data were run on both the FL027 program and the Douglas 
potential flow program [Ref. 1]. The data generated by both 
programs was compared to wind tunnel data for the NACA 572 
wing [Ref. 81. The results are presented in Fig. 5.1 as 
plots of lift coefficient versus a ngle-of-a ttack . The 
results shew that for the NACA 572 wing the FL027 program 
mors accurately predicts the wing lift coefficient than does 
the Ecuglas nregram. 

C. AE-4501 CLASS PROJECT 

The final test phase was conducted by introducing the 
FL027 program into the AE-4501 course as a class project. 
This was accomplished to determine student problems/comments 
concerning the data input program FL027IN and to test an 
additional wing shape. The wing chosen for study was that 
of the A-7 airplane. The A-7 wing has a distinct leading 
edge notch at the approximate mid-span. When the planfcrm 
geometry was run with the notch included the FL027 program 
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Figure 5.1. Prograa Calculated and Wind Tunnel Data 

ra n tc completion but gave negative values for section and 
total induced drag coefficient. The value for the lift 
coefficient was low for the freestream Mach and 
angle-of-attack used. It was found that if the notch was 
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excluded from the wing geometry input data the program 
results were satisfactory both for induced drag and lift 
coefficient. 

From the AE-4501 class experience it was determined that 
sharp wing planform discontinuities cannot be handled by the 
program. If however, the changes in shape are gradual, such 
as a wing glove, the program output appears to be 
satisfactory. Such was the case with the acceptance test 
case data where the wing geometry was that of the F-8 



supercritical wing which incorporates a wing glove. 
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C This JCL routine allocates sufficient space on the mass 
C stcraqe svstem to store the entire tape contents 
//JACK JOE (3266,0178) ,' PASCHAL L-275 9 • ,CLASS=A 
//* MAIN ORG=NFG VM1 . 3 266 P 
// EXEC PGM=IEFBR14 

//DD 1 DD UNIT = 3 3 30 V , MS VG P=PU B4C,DISP=(NEW f CATLG) , 

// DSN=MSS. S3266. WF LOW . DATA, S P ACE= (CYL , ( 16,4,2)) 

/* 

// 



C This JCL routine is used to t 
C a partitioned data set in the 
//JACK JOB (3266,0 178) ,' PASCHAL 
//♦MAIN OEG=NPGVM1. 3266P 
//COPY P HOC FI1E=,MEM= 

// EXEC PGM=IEBGENER 
//SySPRINT DD S YSOUT = A 
//SYSIN DC DOMMY 
//SYSUT1 DD 0NIT=3 4 0 0-6 . V CL=S EE 
// LAB EL= (&FILE ,B LP , , IN ) , 

// DCB= (RECFM=F,BLKSIZE=80, DEN 
//SYS0T2 DD DISP=(OLD, KEEP) , 



ransfer all 
mass stora* 
L-2759 ' ,CLA: 



tape files 
f e system 
;s=j 



to 



// 



DSN = M SS . S3266. WFLOW. SOURCE ( 



// 


DCB= (RECFM ; 


// 


PEND 




// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 


// 


EXEC 


COPY 



=WFLOW, DIS P= (OLD, PASS) , 

=3 , OPTCD=Q ) 

&MEM) , 

IZE=64 00) 



COPY , FIL E= 1, MEM=FL027 
COPY, FILE=2,MEM=A4 111 N 
COP Y , FILE= 3, M E M=VW IN 
COPY,FILE = 4,ME M=A4 1 1 A 01 
COPY , FILE=5, MEH=INTER P 
COPY, FILE = 6, ME M=A4 1 1 P S 
COPY , FILE= 7, MEM=A4 IIP 1 
COPY , FILE = 8, MEM=A4 1 1 P 2 
C0PY,FILE=9. MEM=COUPLE 
COPY , FILE = 10, MEM=ITER 



COPY , FILE= 11 , M E M= D AT A 
COPY , FILE= 12 , MEM=FINA 
COPY, FILE= 13, M EM=BOEB 
COPY, FILE=14,MEM=CONT 
COPY , FILE= 15, MEM=CORD 
COPY, FILE = 16 , MEM=S THE 
//CCPY2 PROC FI LE= , MEM= , LEECL=8 
// EXEC ?GM=IEBGENER 
//SYSPRINT DD SYSOUT=A 
//SYSIN DD DUMMY 

//SYSUT1 DD UNIT =3 40 0-6, VOL=S ER = WFLOW, DIS P= (OLD, P ASS) 
// LABEL= (&FILE,BLP ,, IN) , 

// DCB= (RECFM=F,BLKSIZE=6LRECL 
//SYSUT2 DD DIS?=(OLD, KEEP) , DSN 



IN 

L 

1 

' PLT 
• PLT 
PLT 

0, ELK=6400 



// 

// 

// 

// 

// 

// 

/* 

// 



DCB= (RECFM=FB, LRECL = SLRECL, 
PEND 

EXEC COPY2,FILE=17,LRECL=1 
EXEC COPY2 ,FI LE=18 , LEECL=1 
EXEC COPY2,FILE=1 9, LRECL=1 
EXEC COPY2,FILE=20, LSECL=1 



, DEN=3 ,OPTCD=Q) 
=MSS.S3266.WFLOW.DATA (SMEM) , 
BLKSIZE=&BLK) 

50,BLK=6000,MEM=OUTF27 
50, ELK=600 0, M EM=OUTIFC 
50 , BLK=60 0 0 , MEM=0UT4 1 1 L 
50, BLK=6000,MEM=OUT41 1 U 
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from mass 
accessed by 



C This JCL routine moves all source code files 

C storage to the MVS 004 disk which can be ^ 

C entering GET MVS then following the screen instructions 
C to move source files to your disk. If you want to move 
C the data files to MVS 004 then change the word SOURCE 
C to DATA in the JCL program below. 

//JACK JCB (3266,0178f • PASCHAL L-2759’ , CL ASS=A 
//♦MAIN OR G=NPGVM1 . 3 266 P 
// EXEC FGM=IEBCOPY 
//SYSPRINT DD SYSOUT=A 

//FROM DD DISP = SHR,DSN = HSS.S3266.WFLOW. SOURCE 



//INTO DD UNIT=3350 , VOL=SER=MVS 004, DISP= (NEW, KEEP) , 
// SP ACE= (CYL, {16, 4, 10) , RLS E) ,DSN=S32 66. SOURCE 

//SYSUT3 DD UNIT=SYSDA, SEACE=(CYL, (2,2) ) 

//SYSUT4 DD UNIT=S YSDA, S EACE= (C YL, (2,2)) 

//SYSIN DD * 

COPY OUTDD=INTC, I NDD=FPCM 

/* 

// 
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APPENDI X B 



LIEBARY OF AIRFOIL SECTION GEOMETRIES 

= user input section coordinate data 
= flat plate data 

= symmetrical wing (115? thickness at 30% chord) 

= supercritical wing (cambered, 12% thickness at 32% 
chord) 

= NACA 24-30-0 (cambered, 12% thickness at 30% chord) 

= F-14 wing (cambered, 9.5% thickness at 37% chord) 

= A-7 wing (7 deg droop at 20% chord, 7% thickness at 
43% chord) 

= LISSAMAN 7769 Airfoil (cambered, 11% thickness at 30% 
chord) 

= NACA 0010 (symmetrical, 10% thickness at 30% chord) 



= 


NACA 


00 10-34 


(symm et rica 1, 


10% 


thickness at 


40% 


chord) 


= 


NACA 


0010-35 


(sy mmet rical. 


10% 


thickness an 


5 0% 


chord) 


= 


NACA 


0010-64 


(sy mmet rica 1, 


10% 


thickness at 


40% 


chord) 


= 


NACA 


0010-66 


(symmetrical. 


10% 


thickness at 


6 0% 


c hord) 


= 


NACA 


16-009 


(symmetrical , 


9% t 


hickness at 50% chord) 


= 


NACA 


63-010 


(symmetrical. 


1 0% 


thickness at 


35% 


chord) 


= 


NACA 


63 A0 10 


(symmetrical. 


1 0% 


thickness at 


35% 


chord) 


= 


NACA 


64-010 


(symmetrical. 


1 0% 


thickness at 


40% 


chord) 


= 


NACA 


64A010 


(s ym metr ical , 


1 0% 


thickness at 


40% 


chord ) 


= 


NACA 


65-010 


(symmetrical. 


10% 


thickness at 


40% 


chord) 


= 


NACA 


65A010 


(symmetrical , 


1 0% 


thickness at 


40% 


chord) 


= 


NACA 


66-010 


(symmetrical , 


1 0% 


thickness at 


45% 


chord) 
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